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Abstract 

The non-trivial electronic topology of a topological insulator is so far known to display signatures 
in a robust metallic state at the surface. Here, we establish vibrational anomalies in Raman spectra 
of the hulk that signify changes in electronic topology: an phonon softens unusually and its line- 
width exhibits an asymmetric peak at the pressure induced electronic topological transition (ETT) 
in Sb2Se3 crystal. Our first-principles calculations confirm the electronic transition from band to 
topological insulating state with reversal of parity of electronic bands passing through a metallic 
state at the ETT, but do not capture the phonon anomalies which involve breakdown of adiabatic 
approximation due to strongly coupled dynamics of phonons and electrons. Treating this within 
a four-band model of topological insulators, we elucidate how non-adiabatic renormalization of 
phonons constitutes readily measurable hulk signatures of an ETT, which will facilitate efforts to 
develop topological insulators by modifying a band insulator. 
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Topological insulators (TIs) are a new class of materials which exhibit an electronic 
band gap and a topologically nontrivial electronic structure in their bulk form [l-5], of 
which the latter has interesting consequences for electronic states on their surfaces; for 
example, one-dimensional edge spin states in 2-dimensional quantum spin Hall systems. A 
3-dimensional TI with time reversal symmetry {2!, 
the transport properties of a TI are theoretically predicted to be influenced by the topology 
of its electronic structure jisl, much of its experimental confirmation comes from studies 



6l-ll2| exhibits gapless surface states. While 



of the surface electronic structure [2|], specifically on its odd number of Dirac cones. The 
nontrivial topology of electronic structure results in an intrinsic magneto-electric coupling 
[1^ . yet to be explored experimentally. Although the signatures of the electronic topology 
in the bulk properties may be subtle, they are expected to be more readily detectable when 
there is a sharp change in the electronic topology, i.e. at an electronic topological transition. 

Electronic structure of a TI is characterized by topological invariants s], [l^ , which are 
determined from geometric properties of electronic states as a function of the Bloch vector. 
Similar ideas are involved in the theory of electric polarization 16|], which is determined as 
a geometric phase of Bloch states. An integer quantum of change in polarization has been 
shown to arise from a cyclic evolution (adiabatic pumping) of the insulator along a path 
that encloses a metallic state Similarly, a bulk state with vanishing electronic gap 

in the vicinity of a TI influences the geometric properties of electronic states and is relevant 
to an electronic topological transition (ETT) [l9|. It will be interesting to probe effects of 
the resulting slow dynamics of electrons at the ETT, particularly on phonons which would 
require going beyond the adiabatic approximation. 

An occurrence of broken adiabaticity has been seen in graphene, also characterized by 
a vanishing band-gap: explanation of vibrational signatures of doping in graphene probed 
by Raman spectroscopy required going beyond the Born-Oppenheimer (adiabatic) approxi- 



mation 
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2l| . Similar nonadiabatic effects arising from the vicinity of a TI to a metallic 



state prompt us to carefully look for subtle reflections of an ETT in Raman spectra. 

Among the 3-dimensional time-reversal invariant topological insulators, group V-selenides 



and tellurides (Bi2Te3, Bi2Se3 and Sb2Te3) have been extensively studied |6l. l22|-|27j . In this 
family of compounds, the spin-orbit coupling (SOC) competes with crystal-field splitting 
and chemical bonding (hybridization) between cation and chalcogen, and is strong enough to 
give rise to a topological insulator characterized by an overall odd parity of the valence band 
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manifold (6|]. Vilaplana et al 25|] examined an ETT in Bi2Se3 as a function of pressure using 
Raman spectroscopy and calculations, though the nature of changes in electronic structure 
or topology needs further clarification. Interestingly, a close cousin of these compounds, 
Sb2Se3 is a band insulator with an overall even parity of the valence band manifold, in 
contrast to an overall odd parity exhibited by the TIs (e.g. Bi2Te3) 6|. Under uniaxial 
strain, Sb2Se3 has been shown theoretically to undergo a transition to a TI state {28]. Thus, 
Sb2Se3 is an ideal host system for exploration of a possible ETT and its reflections in bulk 
properties. Competition between the spin-orbit coupling and hybridization can be tuned to 
introduce an ETT either by substitution of Bi at the Sb site or by application of pressure. 
While the former would alter the SOC, the latter would tune the crystal field splitting and 
hybridization between Sb and Se. 

Here, we use a combination of Raman experiments on single crystal Sb2Se3 and first- 
principles density functional theory (DFT) of electronic structure to establish the presence 
of an ETT in Sb2Se3 as a function of pressure. As these calculations do not capture the 
observed phonon anomalies at the ETT, we derive symmetry invariant forms of electron- 
phonon coupling, and go beyond the adiabatic approximation using a model Hamiltonian 
in our analysis, and uncover mechanisms responsible for the anomalies in Raman spectra of 
Sb2Se3 signifying the change in electronic topology at its ETT. In particular, the asymmetry 
in observed phonon anomalies at the ETT arises from reversal of parities of electronic states 
across the gap. 

High pressure Raman experiments on single crystal of Sb2Se3 were carried out at 
room temperature (details in 29|]). Group theory for a centrosymmetric rhombohe- 
dral structure of Sb2Se3 (-Df^) predicts 12 optical zone center phonons represented by 
2Aig (R)+2Eg(R) +2 A2m(IR)+2Eu(1R) symmetry, where R and IR refer to Raman and in- 



frared active modes respectively 25|, |30|. Recently Raman studies have been pursued at 



high pressure to probe the reconstruction of Fermi surface topology giving rise to an ETT in 

-|26|. Probing of ETT with pressure has been done either 



TIs Bi2Te3, Bi2Se3 and Sb2Te3 [22 

by observing changes in pressure coefficients (^), in full width at half maximum (FWHM) 

2^, a-Bi2Se3 
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of Raman bands, or appearance of new modes, as observed in Q;-Bi2Te3 
I and Q;-Sb2Te3 



26j. It is clear from all the above mentioned studies that ETT occurs for 



all these TIs in ~ 3 to 5 GPa hydrostatic pressure regime. Here, we report high pressure 
Raman studies of single crystal Sb2Se3 upto 24.6 GPa to show a ETT transition at ~ 2.5 
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GPa marked by (i) a large softening of the low frequency Raman mode by ~ 16% together 
with an anomalous increase of its linewidth by ~ 200% within a narrow pressure range of 
to 2.5 GPa and relatively less softening by other two modes (M2 and M3), (ii) change in ^ 
and (iii) an appearance of a new mode. 

Figure [1] shows Raman spectra at a few pressures from 0.7 to 20.0 GPa. Three modes 
marked Ml (~ 166 cm~^), M2 (~ 192 cm~^) and M3 (~ 214 cm~^) are clearly seen at 
0.7 GPa. As we shall see later, the mode Ml is assigned as E^, M2 as A^^ and M3 as 
a combination mode. Mode frequencies extracted from Lorentzian fits to the spectra are 
plotted versus pressure in Fig. [2|^a). The data are fitted to the linear equation (sohd line) 
Up = coq + (^)-P in different high pressure regimes marked by vertical dotted lines in 
Fig. [2](a). The data in Fig. 2(a) have been shown to be reproducible in different runs. 

The following observations from Fig. |2]J^a) and (b) are noteworthy: (i) Ml phonon mode 
shows a very large softening of ~ 25 cm~^ upto 2.5 GPa and then the sign of the slope 
(^) changes; whereas M2 and M3 show smaller softening but with a similar change of 
slope, (ii) The full width at half maximum (FWHM) of the Ml mode is maximum near 2.5 
GPa and decreases on either side asymmetrically (Fig. 2(b)), (iii) A new mode M4 appears 
at 2.5 GPa, (iv) at ~ 12 GPa, both M2 and M3 show a slope change, (v) a new mode 
M5 appears around 20 GPa with frequency increasing linearly upto the observed maximum 
pressure (24.6 GPa). The emphasis of the present work is only on the ETT transition at 
~ 2.5 GPa. A transition at low pressure around 3 GPa has been indicated in isostructural 



compounds Bi2Teg 
and their FWHM 



, Bi?Se.s and Sb2Te3 by change in pressure coefficients of the frequencies 
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26[. In comparison, a large softening of phonons upto 2.5 GPa and 
anomalous linewidth maximum of the Ml mode as seen in this work is unique to Sb2Se3. 
Further, we clearly observe three phase transitions as shown in Fig. [2](a). Comparing our 
results with the previous X-ray and Raman studies under pressure on Bi2Te3 [2Ji|], Bi2Se3 



25 | and Sb2Te3 |26|, the transition at ~ 12 GPa can be associated with the rhombohedral 
to monoclinic C2/m phase and transition at ~ 20 GPa to monoclinic C2/c structure. High 
pressure X-ray studies are required to confirm the high pressure phases reported here. 

We now present results of first-principles DFT-based calculations of electronic structure 
and phonons 29|. For the rhombohedral structure of Sb2Se3 with D^^ point group, we 
keep c/a ratio fixed and examine properties as a function of pressure by varying the lattice 
constant a (Figs. HJ^c) and 3). Our theoretical estimate of the lattice constant is a = 4. 09 A, 
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at which a direct energy gap at F-point separates valence band states with odd parity from 
the conduction band states with even parity, as was shown by Zhang et al for a band 
insulator jo]. The odd parity of the state at the top of the valence band is evident in 
the involvement of the p— orbitals of Se atoms present at the inversion center. A strong 
hybridization between p-states of Sb and Se is evident in the electronic eigenstates across 
the gap (Fig. [3]). Band gap vanishes precisely at Pc=2 GPa ( Weyl metallic state) and opens 
up for P > Pc with energies varying linearly with P near Pc. This involves a reversal of the 
parity of valence and conduction bands across Pc marking a transition from a band insulator 

n 

to a topological insulator [19|. Highly accurate calculations using two different techniques 
within adiabatic density functional theory - linear response and frozen phonon methods - 
yield phonon frequencies within a few cm~^ of each other. However, calculated frequencies 
of all phonons at the F point vary linearly with pressure, and exhibit no anomaly at the 
transition at Pc (Fig. S5 in 29|), leaving us puzzled about the mechanism of the sharp 
Raman anomalies. 

To this end, we now present theoretical analysis beyond the adiabatic approximation and 
determine dynamical 20|, |21| corrections to phonon frequencies. We use the universal four 



band model developed by Zhang et al [3 



6i written in terms of Dirac matrices and note that 



states at the top of the valence and bottom of the conduction bands are doubly degenerate, 
and their energies are given by Cq ± Mq, where Cq and Mq are the parameters of the model 
Hamiltonian j3|, and Mq < for a topological insulator. Treating Mq = —k{P — Pc) with 
a positive k, we reproduce the electronic structure near the F point close to the transition 
from a band insulator to a topological insulator (Fig. Wic) and 3) as reflected in the reversal 
of bands of opposite parity |29| . 

Electron-phonon coupling is needed to estimate dynamical corrections to vibrational fre- 
quencies, and we now derive their form at the lowest order within the four-band model, 
expressing them in terms of Dirac matrices. Symmetry properties of the Dirac matrices (ir- 
reducible representations (irreps) for the symmetry group of Sb2Se3) have been derived by 
Liu et al jsl. We note that Raman-active and IR-active modes have even and odd parity, all 
of them are invariant under time reversal. An electron-phonon coupling term in the Hamil- 
tonian can be expressed as a product of phonon degree of freedom and a Dirac matrix (or 
its commutator). A symmetry- invariant electron-phonon coupling is obtained by projecting 
this term onto identity representation of the double group of D^^. For the Raman-active 
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Alg mode, which has the full symmetry of the system, the electron-phonon coupling can 
have all the terms in the four-band electronic Hamiltonian multiplied by the mode displace- 
ment UAig- We pick the leading term that gives changes in band structure obtained by 
distorting the structure with this mode (see inset of Fig. H^): H^^g = A^-^^ua-^^T^, where 
Fs is a Dirac matrix. Similarly, coupling of an IR-active mode A^^ with electrons takes a 
form: Haz^ = ^A2^'UA2u^ i^i where Fjj is a Dirac matrix commutator [Fj,Fj]/2i. Finally, 
symmetry consideration of the lowest order coupling of the E"^ mode with electrons requires 
one to have a quadratic form of Dirac matrices: = AEgiu^^g^i^ + '^\g^2^)^z5i where x 
and y denote cartesian components of displacements of the doubly degenerate E"^ mode. In 
this term, each of matrix F35 and the expression u%Ti^ + m|;_^F25 transform according to 
F^ irrep and is even under time reversal, and their product is symmetry- invariant. Since 
the product of two non-commuting Hermitian operators is not Hermitian, we construct a 
symmetry-invariant Hermitian form: H^Eg = A^;^ [(m|;^Fi5 + u\ F25)F35 — /i.c.]/2i, where h.c. 
means Hermitian conjugate 3ll |. 

We validate the form of electron-phonon coupling by comparing the electronic structure 
of the model with that obtained from first-principles for the lattice distorted with each mode 
of Sb2Se3 (insets of Figs. and b for the case of A\^ and A\g modes). For a frozen A\g mode 
at P = Pc; band gap opens up in a way similar to how it opens up at infinitesimally small 
deviation in pressure or strain. For a frozen A\^ mode, band splitting is more interesting: 
each of the doubly degenerate valence and conduction bands of definite parity split up as 
this mode breaks the parity symmetry. The minima or maxima of bands shift away from 
the F point. However, the splitting or changes in the model band structure associated 
with the electron-phonon coupling of E'^ {^Eg) niode are not captured within first-principles 
calculations, suggesting that the physics of this coupling is beyond the formal mean-field 
description of density functional theory used here. Thus, in the analysis below, we use the 
symmetry invariant form with A-coupling as a free parameter as the other terms are not 
needed to capture the essential aspects of the observed phonon anomalies. 

We obtain the dynamical corrections to phonon frequencies as a function of pressure 
using these forms of electron-phonon couplings (keeping strength of the coupling as free 



parameters) in first-order time dependent perturbation analysis 20|, l29| of the four-band 
model of Bi2Se3 derived by Zhang et al ^ with pressure dependent Mq (see Fig. H]). To 
simplify our analysis, we make use of the layered nature of Sb2Se3 and carry out Brillouin 
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zone integrations only in the ab-plane. It is evident that dynamical corrections to frequencies 
of Raman active modes change sharply below Pc and asymmetrically around Pc. Indeed, the 
sharp drop (Fig. H]) in frequencies of Eg and Afg modes just below Pc, and a gradual increase 
in their frequencies for P > P^. are consistent with the pressure-dependent behavior of modes 
Ml, M2 and M3 in Raman spectra reported here. In contrast, dynamical corrections to IR- 
active mode exhibit a sharp jump above Pc\ Noting the calculated frequencies, we assign 
mode Ml to E| and mode M2 to Alg irreps. The Mode M3 is the highest in frequency, 
and is most likely associated with the sum of two E"^ modes (a second order Raman mode). 
Finally, the new mode M4 appearing above Pc may correspond to (148 cm^^), which 
becomes Raman active due to a change in the overall parity of the occupied states from even 
to odd. 

The linewidth of E"^ mode estimated from our analysis (Fig. 4d and Eq. 9 of [29]) peaks 
asymmetrically near the transition pressure, quite consistent with the observed linewidth 
anomaly seen in our experiments (see Fig. 2(b)). This further corroborates our theoretical 
analysis, and allows us to determine the origin of asymmetry in phonon anomalies to the 
parity reversal of occupied and unoccupied bands in immediate vicinity of the ETT. The 
similar analysis predicts the asymmetry in anomaly of IR active mode contrasting that in 
the Raman mode-with a jump in frequency above Pc- 

In conclusion, we reveal a pressure-induced electronic topological transition (ETT) in 
single crystal of Sb2Se3, a band insulator. A combination of experiment, first-principles cal- 
culations and theoretical model-based analysis presented here show a breakdown of adiabatic 
approximation at the ETT. We established that electron-phonon coupling of nontrivial forms 
leads to anomalies in the observed Raman spectra. These ideas are applicable to electronic 
transitions in other topological insulators too, and expected to stimulate experiments for 
exploring anomalies in IR vibrational spectra, and guide materials scientists in transforming 
a normal insulating material to a topological insulator. 

AKS acknowledges the funding from Department of Science and Technology, India. 
UVW acknowledges funding from DAE-SRC Outstanding Researcher grant. PG thanks 
funding from AFOSR-MURI and an NSF Career award, DMR-0449969. AB thanks CSIR 
for a research fellowship. We thank Prof. M. Cardona for critical comments. 



7 



[1] J. E. Moore, Nature(London) 464, 194 (2010). 

[2] M. Z. Hasan, and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). 

[3] C.-X. Liu, X.-L. Qi, H. Zhang, X. Dai, Z. Fang, and S.-C. Zhang, Phys. Rev. B 82, 045122 
(2010). 

[4] C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 95, 146802 (2005). 

[5] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 106803 (2007). 

[6] H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, Nature Phys. 5, 438-442 
(2009). 

[7] L. Fu, and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 

[8] X.-L. Qi, and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011). 

[9] S. R. Park, W. S. Jung, C. Kim, D. J. Song, C. Kim, S. Kimura, K. D. Lee, and N. Hur, Phys. 
Rev. B 81, 041405(R) (2010). 
[10] S. Mao, A. Yamakage, and Y. Kuramoto, Phys. Rev. B 84, 115413 (2011). 
[11] Y. Xia et al.. Nature Physics 5, 398 (2009). 

[12] W. Zhang, R. Yu, H.-J. Zhang, X. Dai, and Z. Fang, New J. Phys. 12, 065013 (2010). 
[13] D. Culcer, Physica E 44, 860-884 (2012). 

[14] A. Essin, A. M. Turner, J.E. Moore, and D. Vanderbilt, Phys. Rev. B 81, 205104 (2010). 

[15] J. E. Moore, and L. Balents, Phys. Rev. B 75, 121306(R) (2007). 

[16] R. D. King-Smith, and D. Vanderbilt, Phys. Rev. B 47, 1651(R) (1993). 

[17] U. V. Waghmare, Ph.D. thesis, Yale University, 1996. 

[18] K. M. Rabe, and U. V. Waghmare, Ferroelectrics 136, 147 (1992) 

[19] M. S. Bahramy et al. Nature Communications 3, 1159 (2012). 

[20] S. Pisana, M. Lazzeri, C. Casiraghi, K. S. Novoselov, A. K. Geim, A. C. Ferrari, and F. Mauri, 

Nature Materials 6, 198 (2007). 
[21] A. Das et al.. Nature Nanotechnoogy 3, 210 (2008). 

[22] A. Polian, M. Gauthier, S. M. Souza, D. M. Triches, J. C. de Lima, and T. A. Grandi, Phys. 

Rev. B 83, 113106 (2011). 
[23] R. Vilaplana et al, Phys. Rev. B 84, 104112 (2011). 

[24] G. K. Pradhan, A. Bera, P. Kumar, D. V. S. Muthu, and A. K. Sood, Solid State Commun. 



8 



152, 284-287 (2012). 
[25] R. Vilaplana et al., Phys. Rev. B 84, 184110 (2011). 

[26] O. Gomis, R. Vilaplana, F. J. Manjon, P. R.-Hernandez, E. P.-Gonzalez, A. Munoz, V. Kucek 

and C. Drasar, Phys. Rev. B 84, 174305 (2011). 
[27] C. Zhang, L. Sun, Z. Chen, X. Zhou, Q. Wu, W. Yi, J. Guo, X. Dong, and Z. Zhao, Phys. 

Rev. B 83, 140504(R) (2011). 
[28] W. Liu, X. peng, C. Tang, L. Sun, K. Zhang, and J. Zhong, Phys. Rev. B 84, 245105 (2011). 
[29] See Supplemental Material at [ ] for details. 

[30] The nominal frequencies for these modes, for example in Bi2Se3 are Aj^: 74 cm~^, A^^: 174 

cm-\ E^: 42 cm-^ E^: 131 cm"!. 
[31] Similarly, there are other forms of electron and Eg phonon couplings possible: BEg[{u^^r24 — 

</i4)r34 - h.c.]/2i + Be^[{u%T2a - ulTu)Ti2 - h.c.]/2i +CE,{k^u%^ - kyu\)T4, where 

the term with coupling B and B' (C) involves products of (F^) representations, each of 

which is odd under time-reversal. 



9 



(a) 



(b) 





1.6 GPa 


M1 


1.4 GPa 


M1 J^\. 


0.7 GPa 
M3 





M5 20.0 GPa 
M3 



140 160 180 200 220 240 

Raman shift (cm"^) 




140 160 180 200 220 240 260 280 300 

Raman shift (cm'^) 



FIG. 1: (color online) -Pressure evolution of Raman spectra. The solid lines are Lorentzian fits to 
the experimental data points. Appearance of new modes are indicated by arrows. 
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FIG. 2: (color online)-(a) Raman shift versus pressure plot. Solid lines are linear fits to the ob- 
served frequencies. The numbers next to the straight lines are the fitted values of ^ in cm^^ /GPa. 
The error bars, if not seen, are less than the size of the symbol. Fig.(b) shows FWHM of the Ml 
mode (solid points) as a function of pressure and the dashed line is drawn as guide to the eyes. 
Fig.(c) shows the First-principles calculations of electronic structure near the gap as a function of 
pressure in the neighborhood of transition (Pc=2 GPa). 
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FIG. 3: (color online)- First-principles evidence for the pressure induced electronic topological 
transition in Sh2Se-i. Iso-surfaces of charge densities associated with electronic states at the top of 
valence and bottom of conduction bands at T point (a) before and (b) after the ETT, demonstrating 
reversal of bands across the transition. Parity of these states is inferred from the specific Se atoms 
contributing to these states and also from visualization of various spinor components of these states 
(see ^IMJ)- 



12 



p<p 



CO 

"E 

i_ 

re 

S 
< 



"c 

s 



mode 




P< P, 



P> P, 



0.0 
"1 -0.3 

3 




; Eg mode 


^ -0.6 






(0 






3-0.9 
^ .1.2- 




■ 



0.5 
0.4 

0.3 
0.2 

0.1 
0.0 



P<P, 



-1 
-M„ (eV) 



P, P > P, 

c c 





Eg mode 


- E 
: 3 




^ / 




■ >5. / 




^ / 
^ / 
5 / 

■ ^ / 


\ (d) 



■M„ (eV) 



-0.4 -0.2 0.0 0.2 
Mo (eV) 



0.4 



FIG. 4: (color online)- Dynamical corrections to frequencies of phonon modes (a) ^fg, (b) 
and (c) Eg as a function of —Mq {i.e. oc (P — Pc))- The corrections are negative and make the 
modes softer near the transition; their asymmetry allows differentiation between the trivial and the 
nontrivial topology of electronic structure on the two sides of the transition. The insets in Fig. 
4(a) and (b) show how electronic structure changes when atomic displacements of a given phonon 
mode are frozen to distort the structure based on first-principles calculations, and are reproduced 
using electron-phonon coupling (see text) within the 4-band model, (d) The calculated linewidth of 
the mode Eg (Ml mode) as a function of Mq i.e. pressure. 
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